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ABSTRACT 


The HH1/2 system has been observed with the two spectrometers on board 
the Infrared Space Observatory. Diffuse [Oil] 158/urn emission indicate the pres¬ 
ence of a Photo-Dissociation Region (PDR) which is only in part due to the 
external Far-UV irradiation from the nearby Orion Nebula. Additional irradia¬ 
tion must be originated locally and we show that the FUV field produced in the 
recombination regions behind the shocks traced by the HH objects is sufficient 
to induce a PDR at the flow cavity walls. 

The analysis of [Ol], [Sill]34.8pm and [Nell]12.8yum lines suggests shock veloc¬ 
ities v s ~100-M40 kms -1 with pre-shock densities between 100 and 1000 cm -3 . 
H 2 pure rotational lines trace a T~600 K gas which is likely to be warmed up 
in slow, planar, molecular shocks. The coexistence in the areas subtended by 
the instrumental beamsizes of shocks at different velocities is consistent with the 
bow-shock morphology of the HH 1 and, to a lesser extent, HH2 as traced by 
optical images. 

Subject headings: Stars: formation - (ISM:) Herbig-Haro objects - ISM: 
individual objects: HH1/2 - ISM: molecules - Infrared: ISM: lines and bands 


USO is an ESA project with the instruments funded by ESA Members States (especially 
the PI countries: France, Germany, the Netherlands and the United Kingdom) and with the 
participation of ISAS and NASA. 
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1. Introduction 

The HH 1/2 is still the proto-typical optical Herbig-Haro (HH) outflow, and has 
been and it is the subject of large number of studies at almost all possible wavelengths, 
sensitivities and angular resolutions (see “Notes” in Reipurth Rcipurth (1999) for nearly 
100 references on this young stellar bipolar flow). Some of the most recent studies, however, 
have concentrated on the molecular gas properties and its relationship to the optical 
emission, particularly in the context of entrainment, shock chemistry, proper motions, 
multiple outflows, multiple ejection events and the ’unified model’ for stellar jets-CO flows. 
It is in this realm that studies of the far-infrared emission are of great importance, since in 
these collisionally excited objects a large fraction of the atomic and molecular gas cooling 
takes place at these wavelengths (see e. g. Saraceno 1998; Molinari et al. 2000). 

A quick reminder of the HH 1/2 bipolar jet properties goes as follows; it lies in Orion 
at ~ 460 pc, with an angular size of ~ 3' on its brightest component but with signatures of 
shock excited emission from a previous ejection event ~ 23' NW from VLA 1, the outflow’s 
embedded Class 0 driving source (Ogura 1995; Pravdo et al. 1985). The proper motions 
indicate flow velocities of ~ 300 — 400 kms -1 for both the atomic/ionic and warm H 2 gas 
(Herbig & Jones 1981; Eislbffel, Mundt & Bohm 1994; Rodriguez et al. 2000; Noriega-Crespo 
et al. 1997), and this sets a dynamical age of ~ 6 x 10 5 yrs for the outflow, but with a major 
ejection event some ~ 7.5 x 10 4 yrs ago, defined by the brightest condensations. Recent 
single-dish radio observations have shown that the J=l-0 CO emission is bounded by the 
ionic/atomic gas (Moro-Martin et al. 1999), thus providing strong support to the idea that 
molecular outflows are driven by the stellar jets themselves (see e.g. Raga 1994). 

The Infrared Space Observatory (ISO, Kessler et al. 1996) observed HH 1/2 with the 
rnid-IR camera (CAM) and the short and long spectrometers (SWS & LWS respectively). A 
partial analysis of circular variable filter (CVF) ISOCAM observations showed the presence 
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of strong ground rotational H 2 emission lines with a morphology quite similar to that of the 
NIR H 2 gas. These data also showed some signatures of strong J-shocks, like [Nell] 12.8/un, 
as expected in high excitation HH outflows (Cernicharo et al. 1999). In this study we 
present the analysis of the ISO SWS and LWS data with an emphasis on the properties of 
shocks in the HH 1/2 outflow and their effects on the immediately surrounding medium. A 
summary of the observations and data reduction techniques is presented in §2, followed by 
a presentation of the results (§3) and their discussion (§4). 


2. Observations and Data Reduction 

The Long Wavelength Spectrometer (LWS, Clegg et al. 1996) was used in its LWS01 
grating mode to acquire full low resolution (R ~ 200) 43-197/jm scans with data collected 
every 1/4 of a resolution element (equivalent to ~0.07/un for A<90/nn, and to ~0.15/mi for 
A A90//m); a total of 25 scans were collected, corresponding to 50s integration time per 
spectral element. 

LWS data processed through Off-Line Processing (OLP) version 7, have been reduced 
using LIA, the LWS Interactive Analysis 2 Version 7.3. The package consists in a set of tools 
to inspect in great detail the data quality at the photocurrent level. We have refined the 
dark currents and gain for each detector; temporal trends in these two parameters were 
also estimated and subtracted/divided. Finally, the data were recalibrated in wavelength, 
bandpass and flux. The absolute flux calibration accuracy for LWS in grating mode is 10% 
(ISO Handbook, IV, 4.3.2). 

The Short Wavelength Spectrometer (SWS, de Graauw et al. 1996) was used in 
the SWS02 grating line profile observing mode, covering 8 lines: Bry, H2 S(l) - S(5), 


2 LIA is available at http://www.ipac.caltech.edu/iso/lws/lia/lia.html 
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[Nell] 12.8/irn and [Sill]34.8/rm, achieving a spectral resolution R~ 900 — 2000 (see Table 
2.2 ISO Handbook 1.1, p 9). The SWS data was processed using OSIA 3 , the SWS 
Interactive Analysis software, that permits to revise the performance of the dark current 
and photometric checks calibration. As in previous work we have used the latest (November 
2000) bandpass calibration hies to produce the final spectra. 

The final steps of data analysis were done for both LWS and SWS data using ISAP, the 
ISO Spectral Analysis Package Version 2.0 4 . Grating scans (for LWS) and detectors spectra 
(for SWS) were averaged using a median clipping algorithm optimized to hag and discard 
outliers mainly clue to transients. Signihcant memory effects were present in LWS detectors 
LW2, 3 and 4, in the form of systematic differences between the two grating scan directions. 
The LWS observations toward all positions were heavily fringed, partly due to the presence 
in the area of the so-called Cohen-Schwartz (C-S) star (Cohen & Schwartz 1979), partly for 
likely extension of the detected FIR emission; standard techniques available under ISAP 
were used to remove these instrumental effects. 


3. Results 

The region observed with the ISO spectrometers is pictured in Fig. 1, where the 
relevant objects are identified and the locations of the LWS and SWS beams are also 
indicated. The C-S star is a ~30 L & TTauri star (Cohen et al. 1984) that was believed 
to be the exciting source of the bipolar how before VLA 1 was discovered by Pravdo et al. 
(1985). 

3 OSIA is available at http://sws.ster.kulcuven.ac.be/osia 

4 ISAP is available at http://www.ipac.caltech.edu/iso/isap/isap.html 
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EDITOR: PLACE FIGURE 1 HERE. 

Line detection was established by visual inspection of the averaged spectra, and 
Gaussian functions were fitted to estimate integrated line fluxes. In order for a line to be 
considered as detected, it had to be consistently present in both scan directions. This is 
particularly critical for lines in the LWS spectral range between ~110 and 160/xrn, where 
detector memory effects were clearly identified (§2). Figure 2 show the lines detected with 
the two spectrometers. 

EDITOR: PLACE FIGURE 2 HERE. 

To estimate uncertainties we multiplied the standard deviation of the fit from the data 
by the width of the spectral resolution element. Computed values are reported in Table 1. 
Line and continuum flux was observed from all LWS pointed positions. The LWS beam 
has a nearly Gaussian shape, so that radiation may be collected from outside the nominal 
FWHM beamsize. At the distance of one beam this contamination will be of the order of 
10-20%, and since the observed positions with the LWS are adjacent, we can in principle 
estimate the true flux coming from within the individual beams. In the case of the HH 1 
and HH 2 pointings, contamination is only due to the flux coming from VLA 1 (the two HH 
objects are too far to reciprocally contaminate each other). Instead, the observed flux from 
VLA 1 will be contaminated from both HH 1 and HH 2. The contamination corrected flux 
for the lines detected with the LWS are reported in italic in Table 1 

EDITOR: PLACE TABLE 1 HERE. 

Giannini, Nisini & Lorenzetti (2001) have recently independently reduced the same 
ISO-LWS data toward HH1/2. Compared to theirs, our reduction includes the full data 
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reprocessing through the LWS Interactive Analysis (§ 2) which is crucial to really have 
a handle on the level of reliability of the calibrated data. By averaging scan directions 
separately we were able to estimate the line fluxes using the more reliable scan direction 
when detector memory effects were identified (§2). Maybe as a result of this higher level of 
analysis, we have not reported as detected in Table 1 a number of lines (mainly CO lines) 
which are instead reported as detected by Giannini, Nisini & Lorenzetti (2001). However, 
for the lines commonly reported as detected the fluxes are in substantial agreement (within 
10%) for the brighter lines; for fainter lines the discrepancy can be as high as a factor two, 
but the relative uncertainties in the flux estimates are also as large. 

Few basic considerations will help set the frame for the subsequent discussion. 
[Cll]158/jm is detected toward all LWS pointings. Since [Cll]158yum line is one of the 
main coolants in FarUV-irradiated gas, it seems natural to postulate its origin in a 
Photo-Dissociation Region (PDR, see e.g. Ticlcns & Hollenbach 1985). Post-shock cooling 
regions may also produce [Cll]158/zm, but this is not relevant in the present case because 
the implied levels of [Ol]63/jm, expected to be at least ten times stronger than [Cll]158yum 
for a wide range of shock velocities and pre-shock densities (e.g., Hollenbach & McKee 
1989), are simply not observed. The origin of [Ol] lines, also detected toward all observed 
positions, is more difficult to assess; [Ol]63yum, definitely the main shock coolant, can be an 
important (and even dominant) cooling line in PDRs too, if the density is n>10 4 cm -3 and 
the FUV irradiation level, expressed in units of 1.6 10 _3 ergss _1 cm -2 , the local interstellar 
radiation held (Habing 1968), is Go>10 2 (Kaufman et al. 1999). There is very little doubt, 
however, that a major contribution to [Ol] emission comes from shocks since [Sill]34.8/zm 
and [Nell]12.8yum, both important coolants in dissociative shocks, are also detected in the 
present observations; a conclusion further supported by the theoretical J-shock modeling 
(Hartigan, Raymond & Hartmann 1987; Hartigan, Morse & Raymond 1994) applied to 
the strong optical [Ol] 6300+6363 A emission along the entire outflow (Solf & Bohm 



1991). So to explain the observed FIR line emission we need to invoke at least two physical 
mechanisms; a PDR for[Cll] 158/irn and shocks for at least part of [Ol]. A detailed analysis 
of the physical conditions of the PDRs in which the [Cll] 158/irn line is generated will 
provide an estimate of the relative fractions of [Ol] lines radiated by the two mechanisms. 


4. Discussion 

4.1. The PDR 

[Cll]158yum line has similar fluxes, within 20%, along the flow; in particular, the line is 
stronger toward HH2 and fainter toward VLA1. On this basis, we exclude both the C-S 
star and the heavily embedded VLAl source as the FUV (between 6eV and 13.6eV) field 
sources responsible for the PDR conditions. Either the FUV held source is external to the 
HH1/2 how, or, if local, must be somehow distributed along the how. We will see below 
that a combination of the two possibilities best explains the observations. 


4-1-1- External irradiation 

The external irradiation hypothesis is supported by the evidence (Mundt & Witt 1983) 
of a significant UV irradiation on the HH1/2 area from the Orion nebula, whose center 
lies ~ 1.5° to the north. This contamination would amount to 30-r50% of the hux levels 
detected toward HH1/2 shortward of ~1500A, decreasing to ~20% at longer wavelengths. 
Based on the IUE measured spectra (Bohm-Vitense et ah 1982; Bohm, Noriega-Crespo & 
Solf 1993) the integrated FUV held clue to irradiation from the Orion Nebula would amount 
to Go ~7, although it should be noted that the data used by Mundt & Witt (1983) do not 
go shortward of ~ 1500A. This value is compatible with the FUV held from the dominant 
Trapezium star Q\ (about Gq ~ 18), with a reasonable amount of extinction along the path 
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to HH1/2. The maximum [Cll] 158/irn line flux in the ISO-LWS beam (assuming complete 
beam filling) that can be produced by this FUV field is ~ 3.5 10” 19 Wcm" 2 for a gas 
density n~ 10 3 cm" 3 (Kaufman et ah 1999). The [Cll] 158/irn flux decreases by about 40% 
for densities a factor 10 higher or smaller, and is almost linear with Go- We conclude that 
about half, at most, of the PDR [Cll]158yum emission measured by ISO-LWS may be due to 
external irradiation by the Orion Nebula. 


4-1-2. Local irradiation 

In a recent paper (Molinari, Noriega-Crespo & Spinoglio 2001) we proposed that 
the [Oil] 158/un flux detected toward the HH objects in the HH 80/81 flow was emitted 
in a PDR situated at the walls of the cavity excavated by the flow, and irradiated by 
collisionally-enhanced 2-photon FLIV continuum emitted by the post-shock ionized regions 
in the HH objects (Dopita, Binette & Schwartz 1982). We propose that this scenario is 
also pertinent to the HH1/2 system. In terms of observables, this model establishes a 
relationship between the free-free radio continuum emitted by the recombination region in 
the HH object, and the [Cll] 158/irn line flux radiated by the PDR at the cavity walls. In the 
optically thin regime, the former can be expressed as (Curiel, Canto and Rodriguez 1987): 


S v = 1.84 • 10" 4 9 2 


1 -o.i 


10 GHz 


t: 


0.45 


n 


O IQ U S 7 


[1 + 3A83v S7 - 2.745] mJy (1) 


9 is the angular diameter of the recombination region and is estimated off-the-plot from 
the 6 cm VLA D-conf. maps of Rodriguez et al. (1990) to be ~10 7/ . T 4 is the electrons 
temperature in units of 10 4 K and we hold it fixed at 1 in these units. v S7 is the shock 
velocity in units of 100 kms" 1 and we fix it to 1 in these units as determined in §4.2.1 
below. n 010 is the pre-shock density in units of 10 cm" 3 and we will keep this as the only 
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free parameter in Eq.(l). 

The [Cll] 158/jrn line flux can be expressed as (Molinari, Noriega-Crespo & Spinoglio 

2001 ): 


Fcii — 


10 7 An 

"~D2“T 


-l2,.QeV 


r HH Vce fc Xc 


j v dv W cm 


-2 


'6eV 


( 2 ) 


D is the distance of the source from the Sun. j„ is the 2-photon process emissivity and 
is a function, among other things, of the square of gas density. f c is the fraction of the 
2-photon continuum emitted by the HH object which is intercepted by the portion of flow 
cavity wall encompassed by the instrumental beam used to measure the [Cll]158/im flux, 
and can be derived with simple geometrical considerations. The radius of the HH object is 
set to 5" (see above). 


A critical parameter in Eq.(2) is the factor r] ce which accounts for the collisional 
enhancement of the hydrogenic 2s level population above the values predicted by pure 
recombination. rj ce can be determined from a comparison of the predicted two-photon 
spectrum with the observed UV continuum from HH objects. Dopita, Binette & Schwartz 
(1982) find values 5.5 < r] ce < 10.6 for the different knots of HH 1, and 2.8 < r] ce < 6.2 
for the different knots of HH 2. Since most of the knots in the two pointings will be 
encompassed by the LWS beam, we will adopt the mean rj ce values. 


Another important parameter is the fraction Xc of the incident FUV field which is 
re-radiated via the [Cll] 158pm line; this parameter is predicted to vary between 0.1% and 
1% (Tielens & Hollenbach 1985). Since the rest of the FUV flux is reprocessed by dust into 
Far-IR continuum, our data should in principle allow us to independently estimate Xc as 
the ratio between the [Cll]158/im line flux and the integrated FIR continuum as revealed 
by the LWS. Indeed, LWS observations clearly reveal significant continuum emission 
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longward of ^50/jrn toward all pointings. The full spectra are shown in Fig. 3, where the 
minispectra from the ten detectors have been stitched together for cosmetic purposes. The 
region is quite complex at these wavelengths and our continuum data do not represent an 
improvement, given the relatively poor spatial resolution, with respect to previous studies. 
Our data are in substantial agreement with Cohen et al. (1984), also considering that our 
pointed positions do not match the emission peaks of the bipolar emission pattern mapped 
at 100/un with the KAO. The peak wavelengths in the SEDs of Fig. 3 imply black-body 
dust temperatures of 25, 32 and 47 K toward HH 2, VLA1 and HH 2 respectively, in 
agreement with the KAO observations of Harvey et al. (1986). The origin of the FIR 
continuum is not easy to assess. Cohen et al. (1984) do not detect 100/un emission from 
HH 2 to a 3a limit of ~12Jy; clearly, a portion of the continuum that we detect toward 
HH2 is due to contamination from VLAl. As concerns HH 1 it is clear form the KAO data 
that a significant portion of the FIR continuum comes from a radio source ^30" S-SW of 
HH 1 and associated with an H 2 O maser (Pravdo et al. 1985). Based on the same dataset, 
it seems also clear that the C-S star does not seem to significantly contribute longward of 
50/un; this is not surprising given its state of relatively evolved PMS star. 

EDITOR: PLACE FIGURE 3 HERE. 

The pointing where a more reliable estimate of the PDR FIR continuum can be derived 
is probably HH 2. Indeed, the HH 1 continuum is severely contaminated by two sources (see 
above), while the VLAl continuum is likely dominated by thermal dust emission in the 
envelope of this Class 0 source. HH2 instead, may only partially be contaminated by the 
VLAl continuum. The [Cll]158/un/FIR ratio toward HH2 yields \ c > 0.005, clearly in the 
range of values expected from the PDR models. 

With this choice of parameters we run our models (Molinari, Noriega-Crespo & 
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Spinoglio 2001) to determine the quantities in Eqs.(l) and (2) as a function of the gas 
density n (in cm -3 ) and shock compression ratio 7^=n/n 0 . The resulting grids are presented 
in Fig. 4 (full-line for HH2, dashed-line for HH1). The figure also shows the location of 
the HH objects based on the 6 cm flux observed by Rodriguez et ah (1990), and the 
[Cll]158/jm observed with the LWS, decreased by an amount ~ 3.5 1CT 19 Wcm -2 which is 
our estimate of the maximum contribution due to external FUV irradiation (§4.1.1). 

EDITOR: PLACE FIGURE 4 HERE. 

The adopted model seems successful in reproducing the observable quantities. Formal 
values of density and compression ratio are [n, 7U] = [(5200,5100), (40,18)] for HH 1 and 2 
respectively. A problem with this scenario, however, is that we seem to produce too much 
FLIV continuum. The predicted emitted flux at the indicative wavelength of 1500A for 
the adopted model parameters is ~ 8 10~ 13 ergss -1 cm -2 A -1 , which should be compared 
with the dereddened observed values of ~ 10 -13 ergss -1 cm -2 A -1 for the two HH objects 
(e.g., Bohm-Vitense et al. 1982). We can tune the model parameters to produce less FUV 
continuum and still fit the [Cll]158/im and 6 cm observations. For example, we may decrease 
the radii of the ionized regions by 50%, which is not unreasonable since we estimated this 
parameter directly from plots of radio maps. Additionally, we can also increase % c ; in low 
irradiation conditions, like in the present case, this parameter can reach values ~ 0.01 
(Ticlens & Hollenbach 1985) which is not in contradiction with our previous estimate as this 
was a lower limit due to the difficulty of assessing the intrinsic FIR continuum of the HH 
objects (see above). With this choice of parameters we still fit the observed [Cll]158yum line 
flux and 6 cm radio continuum, but with densities a factor 2 higher and compression ratios 
a factor 2 lower, producing an emitted 1500A continuum of ~ 3 10 -13 ergss -1 cm -2 A -1 ; the 
discrepancy with UV observations is reduced to a factor 3, which could be partly accounted 
for by the uncertainties in the reddenening corrections applied. 
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As concerns the jet emanating from VLAl, its low ionization (Solf & Bohm 1991) also 
confirmed by the low radio flux measured (e.g., Rodriguez et al. 1990) is not sufficient to 
produce the amount of FUV continuum that would be required by our models to induce 
the [Cll] 158/jrn emission observed toward VLAl. In the absence of significant dust column 
along the line of sight from the HH objects to the proximity of VLA 1, however, the direct 
FUV irradiation from the HH objects onto the cavity walls at the base of the flow would 
be sufficient to justify the [Cll]158/jm line flux observed toward VLAl. The observation of 
a scattered-light component in the [Sll]AA6716,6731 line along the jet (Solf & Bohm 1991), 
and the presence of a biconical nebula seen in scattered optical light at the base of the jet 
(Strom et al. 1985) would indeed seem to strengthen the idea of a relatively dust-free flow 
cavity (see also Davis, Eisldffel & Ray 1994). 


4-1.3. Line cooling 


We are now able to estimate the amount of cooling expected in lines other than 
[Cll]158//m in the PDR conditions discussed above. The external and relatively faint 
FLTV field (G 0 ~ 7) from the Orion Nebula can generate F([Ol]63/jm)~5%F([Cll]158/im), 
at most, for n~ 10 3 cm -3 . Increasing the density will also increase the expected 
[Ol]63/im/[Cll]158/im cooling ratio (Kaufman et al. 1999), but the absolute [Cll]158/im 
flux emitted would also decrease by a similar amount living the absolute [Ol]63yum cooling 
essentially unaltered. As concerns the local FUV field, the Go value for the PDR will be 
dependent on the assumed distance between the HH objects and the flow cavity walls. The 
CO maps of Moro-Martm et al. (1999) reveal a full width of the low-velocity lobes in the 
HH 1/2 outflow of the order of ~ 30", implying Go ~ 40 at the cavity walls and a maximum 
[Ol]63/im PDR cooling ~ 20 4- 30% of the [Cll]158/rm. As noted above, the absolute [Ol] 
cooling is almost independent on the assumed density. 
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Finally, high-J CO and rotational H 2 lines fluxes comparable to the observed values 
can only be generated in energetic (Go > 10 3 ) and very dense (n> 10 6 cm -3 ) PDRs (e.g., 
Burton, Hollenbach & Tielens 1990, which is clearly not the present case. 

4.2. Shock conditions in the HH objects 

4-2.1. Atomic lines 

The shock conditions toward HH 1 and 2 were diagnosed comparing our observations 
(after correcting the [Ol] lines for PDR contamination, see §4.1.3) with the predictions from 
plane parallel atomic shock models. The shock models were calculated using MAPPINGS2, 
a code developed by Binette, Dopita & Tuohy (1985), which has been thoroughly tested 
(Pequignot 1986). A grid of models was created changing the pre-shock density between 
10 2 and 10 4 cm -3 , and the shock velocity between 60 and 140 kms -1 . The rest of the model 
parameters have been fixed to standard values (Molinari, Noriega-Crespo & Spinoglio 2001). 
The [Ol]63/un/[Ol]145/un— [Nell]12.8/mi/[Sill]34.8/un diagram in Fig. 5 shows that shock 
velocities of v s ~100 or ~140 kms -1 , and pre-shock densities between 10 2 and 10 3 cm -3 
would be adequate to describe the observed line ratios. Higher shock velocities seem to 
be ruled out since [Nell]12.8/nu/[Sill]34.8/im ratios largely in excess with respect to the 
observations are produced (e.g., ~10 for v s ~220 kms -1 , using the planar shock models of 
Hartigan, Raymond & Hartmann (1987)). We note that the beamsize of the SWS at the 
[Sill]34.8//m wavelength is a factor ~ 1.7 larger than at the [Nell] 12.8/jm wavelength, but 
this would not change the position of the HH objects in Fig. 5 by a significative amount 
even for uniformly extended emission. Noticeably, the estimated pre-shock densities are in 
good agreement with the independent estimates based on the shock-excited PDR model 
(§4.1.2). Those models were computed using v s =100 kms -1 , but we note that the adoption 
of v s =140 kms -1 would simply imply, according to Eq. (1), a right-shift of the model grids 
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in Fig. 4, still fitting the observations (with higher compression ratios, though). 

EDITOR: PLACE FIGURE 5 HERE. 

The comparison with optical line diagnostics is not straightforward given the complex 
morphology of the HH objects which appear structured in several knots, both in optical 
and near-IR images, which our observations do not spatially resolve (Herbig & Jones 1981; 
Noriega-Crespo & Garnavich 1994; Davis, Eisloffel & Ray 1994). The SWS pointing toward 
HH1 encompasses knots A, C, D and F; optical lines suggest shock velocities ~ 80km s -1 
for A, and consistently above 130km s -1 for all the others (Hartigan, Raymond & Hartmann 
1987). The dust scattered emission requires at least 150 knis -1 shocks to explain the optical 
observations (Noriega-Crespo, Calvet & Bohm 1991). Our estimate of v s ~ 100km s -1 
might represent average physical conditions over the SWS beamsize, although the negative 
detection of [Olll]88yum in the present observations (expected to be at the same flux level of 
the [Ol]63/zm line) cannot in principle be explained with v s >130kms -1 . The presence of 
strong [O III] 4959+5007 A, however, in HH IF and HH 2A+H (Solf & Bohm 1991) sets 
a lower limit at ~ 100 knis -1 for the shock velocity and suggests that other mechanisms 
like, e.g., partial ionization of the pre-shock material, may be responsible for the depression 
of [Olll]88/jm flux levels. Likewise, the SWS pointing toward HH2 includes knots A, B, 

C, G, H and M, where shock velocities estimated from optical lines range from 100 to 200 
kms -1 (Hartigan, Raymond & Hartmann 1987; Solf & Bohm 1991). 


4-2.2. Molecular lines 

Pure rotational H 2 lines in the ground vibrational level are normally analyzed under 
LTE conditions (e.g., Gredel 1994). Einstein coefficients and wavenumbers were taken 
from Black & Dalgarno (1976) and Dabrowski (1984), and an ortho/para=3 has been 
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assumed. The SWS aperture at the S(l) and S(2) line wavelengths is a factor ~ 1.7 larger 
compared with the other H 2 lines; fluxes for these two lines were corrected accordingly, 
assuming extended and uniform emission. We dereddened the H 2 line fluxes using the visual 
extinctions given by Bohm et al. (1987) and using the Rieke & Lebofsky (1985) extinction 
curve. The H 2 gas temperature and column density can then be directly obtained from the 
slope and intercept of a linear fit to the Boltzmann plots shown in Fig. 6, and are reported 
in Table 2 together with lcr uncertainties in parenthesis; a solid angle equal to the SWS 
focal plane aperture for the S(3)-S(5) lines was assumed, i.e. 14" x 20" or 6.6xlO -9 sr (valid 
for all detected H 2 lines but S(1), where the solid angle is ~30% higher), with a beam Filing 
factor of 1. The fits show that a single temperature component is adequate to represent the 
observations, given the observational uncertainties. 

EDITOR: PLACE FIGURE 6 HERE. 

LIsing the same model and the parameters determined from the Boltzmann plots, we 
can estimate an H 2 cooling, for the component traced by the (0-0) lines, of 9 10~ 3 L 0 and 
2.5 10 -2 L e for HH 1 and HH2, respectively. 

EDITOR: PLACE TABLE 2 HERE. 

The issue of the nature of the H 2 excitation in HH 1 and 2 has been discussed by 
several authors in the last decade using near-IR imaging of the S(l)l-0 2.12/rm and near-IR 
spectra in the H and K bands (Noriega-Crespo & Garnavich 1994; Davis, Eisloffel & Ray 
1994; Gredel 1996; Eisloffel, Smith & Davis 1996; Davis, Smith & Eisloffel 2000). The 
spectra suggest excitation temperatures ranging between 2000 and 3000 K in the various 
emission knots revealed in the images. Detailed modeling suggests that C-type (Draine 
1980) shocks best fit the observational scenario; H 2 emission would originate in the tails 
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of the bow-shocks, whose tips are instead bright in [Fell] 1.64pm (Davis, Eisloffel & Ray 
1994), where the shock velocities are low enough for H 2 not to be dissociated. The fact that 
Boltzmann plots based on H 2 ro-vibrational are satisfactorily fitted with single linear fits led 
Gredel (1996) to doubt as to the C shock origin, since a stratification of temperatures up 
to ~ 3000 K should be expected. Eisloffel, Smith & Davis (1996) review the most popular 
shock models and indeed confirm that curved C-type shocks produce a quite shallow trend 
in the Boltzmann plots; besides, the predicted excitation temperatures of the v = 1 lines are 
significantly lower (~ 1000 K) than actually observed. Planar C-type shocks, on the other 
hand, seem to produce a rather abrupt change of slope at upper energy levels ~ 4000 — 5000 
K. Since this is exactly at the separation between the v — 0 lines we detected in the present 
work and the v > 0 lines detected in the near-IR, the combined analysis of mid- and near-IR 
H 2 lines seems to suggest an origin in planar C-type shocks. 

C-type shocks are good emitters of molecular lines in general. Our marginal detection 
of high-J CO lines at flux levels comparable to H 2 (0-0) lines suggests shock velocities 
15 < v s < 20 kms -1 and pre-shock densities no > 10 5 cm -3 (Kaufman & Neufcld 1996). 
However, for the same parameters the model also predict an H 2 O(2 12 -l 0 i) line flux about 
two orders of magnitude higher than actually observed in HH2 (Tab. 1). Unless H 2 0 
vapor is confined in a much smaller area compared to H 2 or CO, there is a clear deficit 
with respect to expectations. This seems unusual for Class-0 objects (Giannini, Nisini & 
Lorenzetti 2001) and resembles more evolved Class-1 systems. The possibility that water 
vapor depletes onto grains in the post-shock cooling region has been verified for HH 7 in 
NGC1333 with the detection of the 62pin ice band (Molinari et al. 1999); we do not detect 
such a feature anywhere in the HH 1/2 system, although the filling factor of the bow-shocks 
present in HH 1 and 2 (e.g., Fig. 1) in the LWS beam is lower than in the case of HH 7 (also 
because of the larger distance). An interesting alternative to explain the water vapor deficit 
could be offered by the marginal detection of the fundamental lines of the OH molecule in 
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HH2. Given the high velocity of the dissociative shocks in HH1/2, significant amounts 
of Far-UV flux can be expected to be irradiated by the post-shock regions (see §4.1); in 
particular, the photons between 1300A and 2000A are capable of dissociating the water 
molecule into H and OH (Andresen, Thissen & Schroeder 2001). Water photodissociation 
by a UV held was also invoked by Spinoglio et al. (2000) to explain the surprisingly high 
OH cooling found toward T Tail. 

Comparing the v s ~20 kms -1 found for the C-type molecular shocks with the v s ~100 
kms -1 deduced for the dissociative J-type shock traced by the atomic lines (§4.2.1), it 
is plausible that the two shock regions are respectively associated with the tails and tips 
of bow-shocks. Such structures are indeed identified in HH 1 and, with a more complex 
morphology, in HH2 (Hester, Stapelfeldt & Scowen 1998). The 2 orders of magnitude 
difference in pre-shock densities implies that significant compression of the pre-shock 
material has taken place in the molecular shocks. This is expected given the very nature of 
the C-type shocks, where the relatively low shock velocity and the presence of transverse 
magnetic field combine to accelerate and pre-compress the upstream material. 1-band 
polarimetry (Strom et al. 1985) indeed confirms the presence of a magnetic field oriented 
along the HH 1/2 flow. 
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Fig. 1.— HST [Sn] image of the HH1/2 region. The crosses mark the pointed positions, 
while dashed circles and full-line rectangles indicate the FWHM beamsizes of the LWS and 
SWS respectively. The position of the C-S star is also indicated. 

Fig. 2.— a. H 2 , [Nell]12.8/im and [Sill]34.8/mi lines detected toward HH 1 (left column) and 
HH 2 (right column). 

Fig. 2.— b. [Ol]63//m lines detected toward the three positions observed with the LWS. 

Fig. 2.— c. Lines detected between 110 and 180/jrn toward the three positions observed 
with the LWS. 

Fig. 3.— Spectral energy distributions observed with LWS toward all pointings. The 10 
LWS detectors have been stitched together to produce a smooth SED. No correction has 
been applied for the reciprocal contamination between adjacent beams. 

Fig. 4.— [Cll] 158/un line flux as a function of 6 cm radio continuum flux expected for 
PDR emission induced by the FUV two-photon continuum from the ionized regions in the 
HH objects. The grids show the loci of points of constant density n (horizontal lines) and 
constant 7?.=n/no compression ratio in case of a recombination region 5" in radius behind 
a shock front in a HH object with v s =100 kms" 1 . The other parameters used are [r] ce , 
Xc]— [(8-0,4.5),(0.005,0.005)] for the dashed and full-line grids (HH 1 and HH2) respectively. 

Fig. 5.— Diagnostic diagram showing the model results for the indicate line ratios; full lines 
are the loci of constant shock velocity (from 60 to 136 kms -1 ), while dashed lines are loci 
of constant pre-shock density (from 10 2 to 10 4 cm -3 ). 

Fig. 6.— H 2 excitation diagrams (Boltzmann plots) for HH 1 and HH2 (empty and full 
circles respectively). Nj/gj is reported as a function of the energy of the level in K. The full 
lines represent linear regressions to the data points. 
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Table 1. Observed Line Fluxes 


Line 

A(/im) 

VLA1 

HH1 

HH2 


SWS Lines 

(O-O)S(l) 

17.03 


0.7(0.3) 

2.8(0.3) 

(0-0)S(2) 

12.28 


2.3(0.5) 

3.1(0.3) 

(0-0)S(3) 

9.66 


3.0(0.6) 

5.8(0.9) 

(0-0)S(4) 

8.03 


1.8(0.8) 

4.7(0.5) 

(0-0)S(5) 

6.91 


3.6(0.9) 

14(2) 

[Nell] 

12.81 


4.7(0.9) 

6.8(0.7) 

[Sill] 

34.82 


5(1) 

6(1) 

LWS Lines 

CO 18-17 

144.8 

2(1) 



CO 17-16 

153.3 



3.0(0.6) 

CO 16-15 

162.8 

6.8(0.7) 

5(2) 

4(2) 



(5.7) 

( 4 ) 

(3) 

OH [3/2,5/2-3/2,3/2] 

119.3 



4(1) 

OH[3/2,l/2-l/2,l/2] 

163.2 



3(2) 

o-H 2 0[2i 2 -1oi] 

179.5 



7(2) 

[Oi] 

63.18 

173(2) 

101(2) 

150(3) 



(lol} 

(86) 

(135) 

[Oi] 

145.5 

5(1) 

4(1) 

4(1) 



( 4 ) 

( 4 ) 

( 4 ) 

[Cn] 

157.7 

79(1) 

90(3) 

82(3) 



(63) 

(Si) 

(76) 









Table 2. H 2 Physical Parameters 


Object 

Th 2 

N(H 2 ) 


(K) 

(10 19 cm -2 ) 

HH1 

670(40) 

1.1(0.4) 

HH 2 

730(20) 

2.2(0.3) 

Note. 

- Based 

on aperture- 


corrected line fluxes (see text). 


















































